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entanglement- can be identified among the eigcnstates. We also discuss the violation of 
Bell inequalities in this system and show the impact of a limited measurement fidelity 
on the detection of entanglement and quantum nonlocality. 
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1. Introduction 

Entanglement is one of the most intriguing consequences of quantum mechanics [I] and 
has long been debated in a case against its completeness [2]. Entanglement manifests 
the nonlocality of quantum mechanics. Bipartite entanglement is well established in 
optical systems, see, e.g., [31 IU E], and has recently been verified in superconducting 
circuits [HI [7]. Quantum information processing uses entanglement as a resource [H [9], 
i.e., the preparation of entangled states is essential for quantum computing. 
In multipartite systems, different types of entanglement can be classified. Specifically, 
in a tripartite system, one can recover the same type of bipartite entanglement as in 
a two-particle system - the measurement of one of the particles reduces entanglement 
by taking out the particle being measured but leaves the entanglement between the 
remaining particles intact. On the contrary, there can also be tripartite entanglement 
where a single one-particle measurement completely destroys entanglement between all 
parties. 

Quantum computing and control is very mature in liquid-state nuclear magnetic 
resonance and in atomic and optical physics. In the former, due to the use of pseudopure 
states, it is not obvious that strong entanglement can be created. Creation of entangled 
states of ions, neutral atoms, and photons [31 01 El El [7] has been highly successful. The 
required experiments, though, become increasingly complex, caused by the very weak 
interaction of the objects to be entangled, which often is also only effective for a very 
short time. 

On the other side, condensed matter systems typically have strong inter-particle 
interactions which can act effectively because the particles occupy fixed positions in 
space, e.g. in a lattice. Thus, entangled states are not that exotic in correlated condensed 
systems [ID] . Going beyond the simplest system, namely the spin singlet, we consider a 
single triangle of spins with antiferromagnetic Ising coupling. Without external field, the 
ground state of the system is spanned by the three degenerate frustrated states in which 
the orientation of one of the spins differs from the other two. Creating a superposition of 
these states right away leads to maximally entangled Werner (W) states [H] which only 
contain bipartite entanglement. These superpositions can be created in the quantum 
version of this model, i.e., if the system is put into a transversal magnetic field. Spin 
Hamiltonians with 'designed' properties can be implemented using macroscopic devices 
such as superconducting qubits rather than elementary spins [121 H3]- In particular, 
superconducting flux qubits allow for a strong inductive qubit-qubit interaction. 
The paper is organized as follows: In section [2] we propose a design for the coupling 
of three flux qubits. It provides a strong interaction whose strength can be designed 
over a large range during fabrication. The form of the eigenstates is discussed in section 
El states in the proximity of Greenberger-Horne-Zeilinger (GHZ) states [H] as well 
as W states can be identified. Section H] describes the characterization and detection 
of tripartite entanglement in this system. In section El we discuss the violation of 
Bell inequalities by virtue of optimally chosen operators and compare the required 
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Figure 1: The design of the flux qubit triangle. The three qubits are formed by the 
three small isosceles triangles, the round arrows in the qubits define the directions of 
the currents. Small crosses represent the Josephson junctions in the individual qubits, 
large crosses the coupling junctions. 



We consider the triangle design sketched in figure [U consisting of three three- Josephson- 
junction (3JJ) flux qubits [15J. The individual qubits operate in the flux regime with 
k B T -C <C Ef 1 <C A, % — 1, 2, 3, and A, Eq, and are the superconducting gap, 
the charging and Josephson energy of the i-th qubit, respectively. Therefore, the two 
different directions of persistent currents in the qubit loop form the basis states of the 
qubits. As this qubit design shows quantum behaviour even for small self- inductance, 
it allows for the fabrication of circuits that are relatively insensitive to external noise. 
However, a small self-inductance results in a very weak inductive coupling. This makes 
simple coupling schemes, like placing disconnected qubits close to each other, insufficient 
as they are restricted by the mutual geometric inductances. In contrast, the coupling 
between the qubits in this design is provided by shared lines with additional large 
Josephson junctions [T6] . 



As the potential energy is dominated by the Josephson energies of the junctions, we 
will only take this energy into account for the calculation of the coupling strength and 
neglect the inductive energy from the geometrical inductance of the triangle and the 
kinetic inductance of the shared lines. For convenience, we will consider the Josephson 
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energy of the qubit junctions (the ones on the outer edges of the triangle) and the large 
shared junctions separately. 

The total Josephson energy of the qubit junctions is given by 

3 

£jos,q = -Ej J^(cos tpf 1 + cos <p$ } + a cos (pf) . (1) 

i=l 

Combining the additional phases across the shared junctions into the phase 0^ (e.g. 
</>W = (p(<p — for qubit 1), the fluxoid quantization reads 

^ + ^ + ^+^ + ^ = 0, (2) 
where $ x denotes the externally applied flux. 

The coupling junctions are large compared to the qubit junctions and their critical 
currents are far above the persistent currents I v in the qubits. Hence, their phases are 
small and we obtain for first order (for qubit 1) 

orW _ r(2) _ A3) 

(i) w h h_ ; (3) 

Ic,§ 

where the critical current Jc,s is assumed to be equal for all three shared junctions. 
In order to separate the effect of (f>^\ we solve the fluxoid quantization for, e.g., tp^\ 

a cos = a cos ( + $ + ^ J + , (4) 

where we again expanded to first order in 0W and used 27r$x /$o ~ vr as well as the 
relation ip^ = ip^ = ±<^* with cosy?* = l/2a [17] for the minima of the potential 
landscape of a single qubit. Equation ([[]) then can be rewritten as 

E, n = Ei /(') 2 + JLf (5) 

-^Jos,q -^Jos,uncp 2e 2 Ei*^ p 2e 2 Ei*^^ p p ' 

J ' b i=i J,b i=i j>i 

with the Josephson energy of the uncoupled qubit system -Ej os ,uncp- This contribution 
to the coupling is thus found to be pairwise antiferromagnetic, with a larger shared 
junction resulting in a weaker coupling. 

The second major contribution to the potential energy is the Josephson energy of the 
shared junctions 0, 

3 3 / (i)2\ 

£jos,S = -Ej, S C0S <Ps ~ "^J,S Z) 1 " • ^ 

i=l i=l \ / 

,(i)~/-r(i) r(2)^ 



Using ~ (/p — Ip )/Ic,s etc-; it reads 

£j „, s = _ 3£j , s t^Ef-^tE « , (7) 

' i=l ' i=l j>i 

X Note that we expand to second order in ip^ here; this is justified by Ej$/Ej = Iq.s/Ic which makes 
© and ((U) of the same order in /c/^c,s- 
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yielding a ferromagnetic coupling with half the strength of the qubit term. 
Hence the total the coupling strength is given by 

AC/W) = C«) = -^/«/W + 0[{E z /E^f]. (8) 

Expressing the persistent currents in the qubits in terms of the Josephson energy of the 
qubit junctions, I p = Iq^I — 1 /4a 2 = 2eEj/h ■ a/1 — l/4a 2 , the coupling is seen to 
depend directly on the size ratio r = Ej/Ej^s of the qubit's and shared junctions, 

In experiments with inductive coupling [16] or shared lines [18] between the qubits, the 
coupling is very small (C/Ej ~ 0.1% — 0.5%), reflecting the small mutual inductance 
(geometric and/or bulk superconducting kinetic inductance) that can be achieved in 
such designs. Contrastly, ([9]) allows for coupling strengths which are one order of 
magnitude stronger. The shared junctions are still large and linear: for a size ratio 
r of e.g. 1% — 5%, the shared junctions basically behave like a linear inductor, 
Aj,s — $0/(2^/0,3 cos 0s) ~ h 2 /(4e 2 £'j j 5). Thus, the shared junctions can be viewed as 
ultra-compact inductors making the qubit loops equivalent to the standard 3-j unction 
flux qubit with a moderate inductance [i~5l [T9] . 

The currents in the qubits are quantum mechanically associated with a z operators and 
the effective Hamiltonian reads in terms of the Pauli spin matrices 

H = E H ei *® - l Ai *®) + E E aW- ( 10 ) 

i=l ^ ' i=l j>i 

The energy bias e can be tuned by the externally applied magnetic flux, whereas the 
tunnel matrix element A depends on the size ratio of the qubits' junctions (TTIED]- In the 
following we assume the qubits to be identical (Ai = A 2 = A 3 = A, e\ = €2 = €3 = e). 
Note that even though the tunnel splittings Aj exponentially depend on fabrication 
parameters [T7], recent experiments [21] have made them extremely close to each other 
and, in principle, their values can be externally adjusted [22]. 

3. Eigenstates of the system 

We aim for preparing tripartite entangled states in a preferably easy and stable way. 
Both demands are naturally met by the eigenstates of a system, as they are easy to 
prepare by 7r-pulse driving and stable against pure dephasing processes. Figure [2] 
displays the eigenenergies as a function of the energy bias e for a given coupling strength. 
These have been found numerically but can be understood analytically. 
For large positive or negative energy bias, the ground state and the highest excited state 
are the classical, totally polarized states |ttt) or respectively, whatever is favored 

by the magnetic field. In this case, the antiferromagnetic coupling can be disregarded. 
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Figure 2: Plot of the eigenenergies of the eigenstates \Ei)-\Es) for the indicated 
coupling strength. \E§) and \E§) as wen as l-^V) an d \E%) are always degenerate. 



The other excited states are accordingly the classical frustrated states with the three 
states | til), | III)) an d lilt) being degenerate, as well as ||tt)> I tit)) an -d |tt!)- 
Going to finite energy bias, we need to take the coupling into account (C = 1.4 A in the 
following). The antiferromagnetic coupling energetically favors frustrated states. For 
e = 0, the ground state \E\) therefore contains a larger contribution of frustrated states 
and a smaller contribution of polarized states, 

\Ek) = J— {\U) + *(lttt> + HID)} ) (11) 

Vo + 2 o z 

where |/ a n) denotes the equal (non- normalized) superposition of all frustrated states 

|/all) = Itl!) + Itlt) + Itt!) + lltt) + lit!) + lilt) , (12) 

and S ~ 0.2 is small, i.e. the aligned states |ttt) an d ||||) a re suppressed. 

The highest excited states \Eg) and \E$) show the opposite behaviour and consist for 

e = mainly of a superposition of distinct polarized states, 

\E 3 ) = — L= {[4^) + |ttt) " ^ll/all)} (13) 

V2 + 6 8? 

with 8\ w 0.07, and 

\E A ) = y J_ {[Ui) - Ittt) + fcl/*)} (14) 

+ 6<5n 2 
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with 62 ~ 0.1. Thus, we find these states to be in the proximity of GHZ states, commonly 
represented by |GHZ) = (\Ui) ± |ttt))/V2. 

Due to the large antiferromagnetic coupling, the eigenstates discussed above evolve 
into classical states already for small detuning of the energy bias. However, there are 
more regimes of tripartite entangled states among the spectrum of eigenstates: at finite 
positive and negative energy bias e = ±e* ~ ±2.6A two more anticrossings involving 
\E 2 ) and \E 3 ) are present. The explicit form of the state forming the lower branch at 
e = — e* reads 

l£ 2 )_ e » = J_ {im> + mt) + litt) - 1111)} + 

+ Ss{\m) + im> + im> - ittt)} , (15) 

where £3 ~ 0.09. Therefore, this state is close to 

|GHZ) = ^(Ittl) + Itlt) + litt) - llll)) = -^=(|000> + |111>) (16) 

with |0) = (|t) + i\i))/V2 and |1) = -(||) - i\l))/V2. As |GHZ) and |GHZ) can be 
transferred onto each other by purely local operations, they have identical entanglement 
properties. 

The dash-dot lines in figure [2] indicate two two-fold degenerate subspaces. It is shown 
in the appendix that arbitrary states in these two-dimensional Hilbert spaces can be 
prepared by coupling these subspaces to the ground state via resonant driving. The 
subspaces are spanned by frustrated states. Since any frustrated state |/) is an eigenstate 
of the coupling operator, (ai^a^ + a^ai^ + a^ai 3 ^)\f) = —\f), the form of an 
eigenstate prepared in these subspaces does not change with the coupling strength. 
Among the eigenstates, W type states can be found, with the maximally entangled 
state (with respect to the measure of global entanglement [23]) in the lower energy 
subspace for e = taking the form 

i^Lx)=^{2(im)+iut))-(i-^)(itit)+im))-(i+^v / 3)(im}+iitt))}.(ir) 

The local transformation U L iV^ax) = |W) rotating iV'max) on t° the common 
representation of a W state, \W) = -^(Ittl) + I tit) + litt))) is given by 

U L = {e^/ 3 Ry(n/2)}^ ® {R z (2ir/3) R y (n/2)}^ ® {R z (-2tt/3) R y (7r/2)}^ , (18) 

where R z {9) (R y (8)) is a rotation around the z-axis (y-axis) by an angle 9. Thus, we 
can see explicitly that the state iV'max) i s indeed locally equivalent to a W state. 
With respect to the occurrence of W states, the relative pairwise coupling strengths 
between the qubits are essential: equal mutual (antiferromagnetic) couplings like in the 
design proposed here make the system frustrated and facilitate W-like eigenstates. This 
is not given in, e.g, a linear chain of qubits where the two nearest-neighbour couplings 
are stronger than the next-nearest neighbour coupling between the outer qubits. 
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4. Detection of tripartite entanglement 

Conjecturing entanglement from the energy spectrum is not a satisfactory experimental 
indication. We thus want to discuss measurement protocols to verify entanglement more 
explicitly. 

A necessary and sufficient condition for the entanglement of a given state is only known 
for two-qubit systems [23]. The determination of entanglement for multipartite states, 
however, is an open question. For three qubits two inequivalent kinds [25J of genuine 
multipartite entanglement (i.e. each party is entangled with each other party) can 
occur, represented by the GHZ state and the W state. A tool for detection of any 
kind of genuine tripartite entanglement for arbitrary states is not at hand; however, if 
some knowledge about the state under investigation is provided, entanglement witnesses 
(EWs) can be used [261 EZ] • These are observables with a positive expectation value for 
all (bi-) separable states (in general n — 1 partite entangled states), whereas a negative 
expectation value indicates the presence of tripartite (n-partite) entanglement. The 
common way to construct an EW W for a state is 

W = al-|^)M , (19) 
where a is the maximal squared overlap of \ip) with any biseparable or fully separable 
state. Determination of a is in general complicated [271 EE1 EH], but we can use the 
aforementioned proximity of the states under investigation to W and GHZ states, 
respectively, to make use of known values for a. We will discuss how to measure EWs 
in the end of this section. 

Figure [3] shows the expectation values for two EWs [E^W-^^E-i) and (^IVV^^-E^) 
for varying energy bias. The EW W^^- is constructed such as to detect states of the 
form given in ( JT5]) . i.e., I-E2) at the left anticrossing. In order to construct the optimal 
EW, one would have to determine the maximal squared overlap of jE^) in this point 
with any non-GHZ entangled state. Instead, we make use of the vicinity of \E2)- e * 
with I GHZ) defined in ([T6j) and use [30] = § 1- | GHZ) (GHZ |. The optimal local 



decomposition of VV^^ requires four experimental settings [31], see |Appendix B[ Since 



the system is invariant under a combined flip of the spins and an inversion of the sign 
of e, VV^^ (^ ne optimal EW at the right anticrossing) is the totally flipped counterpart 

Figure [3] also shows the 3-tangle r [32], which allows for a reliable distinction between the 
two classes of entanglement, as it is zero for all W type states (and all separable states, of 
course), whereas it takes positive values for all states in the GHZ class. However, the 3- 
tangle can only be measured by full state tomography which requires many more settings 
than the witness. Both quantities indicate a strong (limiting case: r max = t(|GHZ)) = 1, 
(W^j^nun = (W^^)min = —1/4) tripartite entanglement of GHZ type over a large 
range of the energy bias. 

In Figure HI the 3-tangle and the expectation value of the corresponding EW for 
l^max)' = § — l^maxX^maJ [32] is displayed. Its expectation value is positive 
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Figure 3: 3-tangle and expectation value of the GHZ witnesses W-=, TI 
the state \E2). For finite energy bias e = ±e* we find a peaking 3-tangle as well as 
negative expectation value for the two GHZ witnesses, indicating entanglement of the 
GHZ type. Moreover, the entanglement is relatively robust to detuning of the energy 
bias. 




-'max) 



Figure 4: 3-tangle and expectation value of the W witness W^y f° r the state 
The vanishing 3-tangle excludes entanglement of the GHZ type, whereas the negative 
expectation value of the W witness indicates a W type entanglement. 
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for biseparable and fully separable states. 

It thus detects genuine tripartite entanglement in general, without distinguishing 
between W and GHZ type entanglement. In connection with the 3-tangle a distinction 
can be achieved, though, stating an entanglement of the W type in a large range for e, 
approaching its theoretical minimum ((W\y) = — 1/3) for zero energy bias. 
Measuring these entanglement witness only requires single-qubit measurements of 
certain linear combinations of Pauli matrices. As a consequence, in a superconducting 
qubit setting, this can be achieved by decoupling the qubits, performing an appropriate 
single-qubit rotation, and then measuring the most accessible variable such as the flux, 
here corresponding to a z . Finding such a local decomposition is a demanding task that 
has been solved in literature [33l [31]. The appropriate optimal local decompositions 



for the witnesses employed here are shown in Appendix B Local measurements of this 
type have been performed in the context of Bell tests 0, [7] . 

We have to note that in the proposed setup, the couplings between qubits are 
permanently switched on. However, this coupling can be effectively removed by 
nonadiabatically pulsing e to strongly detuned values. The evolution in this detuned 
state corresponds to single-qubit rotations alone and thus does not change entanglement. 



5. Violation of Bell inequalities and robustness to limited measurement 
fidelity 

So far, we made use of EWs as tool for the detection of tripartite entanglement. Another 
common approach for the detection of entanglement are Bell inequalities. Multiqubit 
states can contradict local realistic models in a new and stronger way than two-qubit 
states [34], [35], reflected by a stronger violation of Bell inequalities [36], [37]. In the case 
of three qubits, quantum mechanics predicts (GHZ|Mghz|GHZ) = 4 for the expectation 
value of a Bell operator Mqhz in the GHZ state [§|, while the local prediction gives an 
upper threshold of (?/>|Mghz|'?A) < 2 for any states In comparison, the maximal 

expectation value of a Bell operator for a two-qubit state is 2\/2, in contrast to the 
local prediction of < 2, yielding a maximal violation only by a factor of However, 
Bell inequalities are not distinctive to the type of entanglement; moreover, there are 
entangled states that do not violate any Bell inequality. Nevertheless, their violation 
as sign for non-classical correlations is highly substantial as an ingredient to quantum 
information processing. Besides, it also allows for a comparison between the twopartite 
case and the tripartite case with respect to the robustness to limited measurement 
fidelity. 

We again investigate the two states \E 2 ) and |^max) f° r varying energy bias, see figures 
[5] and [6] For \E 2 ), we use the Bell operator [36] Mqez = ®xO~xO~ z + o~ x a z a x + 
&z <?x & x — d- z cr z cr z . As for EWs, we can determine optimal Bell operators for given target 
states. In general, Bell inequalities for three qubits are constructed from the operator 

§ the specific form of Mghz will be given later 
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Figure 5: 3-tangle and expectation value of the Bell operator M GHZ for the state 
l-E^)- A significant violation of the corresponding Bell inequality can be observed. 

rh(a,b,c) = (a ■ a) <S> (b ■ a) (g) (c • a), where a, b, and c are real three-dimensional 
normalized vectors, which define a rotation of the Pauli matrices a = (a x ,a y ,a z ). A 
Bell operator is then given by 

M = m(a, 6, d) + m(a, b', c) + m(a', b, c) - m(a', b', d). (20) 

In order to obtain an optimal Bell operator for a given state, one can optimize over the 
six unit vectors a, a', b, b', c, and d [231 EH]- The optimal values for M w adapted to 
the state |W), as obtained by a numerical optimization, are listed in table [TJ The Bell 

Table 1: Components of the vectors a, a' , b, b', c and d for the Bell operator Mvv- 
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operator for the state |^max) then reads M^v = U L ^ My? U L with £/ L being the local 
propagator from (fTBl) . 

Any experimental test of tripartite entanglement or the violation of Bell inequalities 
involving three qubits will be more fragile than a two particle test and will be 
put in jeopardy by detector imperfections (as three-party correlations need to be 
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e/A 

Figure 6: 3-tangle and expectation value of the Bell operator for the state 
l^max)- The maximal violation of the Bell inequality is not as high as for \E2) in 
figure however, the violation persists over a larger range of e. 



measured in either case, the measurement fidelity enters with the third power) and 
fabrication uncertainties. However, for the Bell inequalities, the stronger violation might 
compensate for that. We consider the effect of a limited measurement fidelity / < 1 on 
the expectation values of the EWs and Bell operators introduced above and compare 
the results to a representative two-particle case. We model a non-perfect measurement 
of a spin component a, by the perfect measurement of a spin component a[ which yields 
the correct measurement result with a probability / and '1' otherwise, 

^ = /a 4 + (l-/)l. (21) 

Figure |7] compares the decay of the Bell violation with decreasing measurement fidelity 
for the cases discussed above and for a representative two-qubit case, namely the 
violation of the Clauser-Horne-Shimony-Holt (CHSH) inequality [39]) by a Bell pair 
|^b) = Tjdtl) ~~ I It})- We find that the larger initial violation for the tripartite 
cases compensates for the quicker decay of the violation and allows for a slightly lower 
minimal detector fidelity; table |2] summarizes the minimal detector fidelities for the 
detection of tripartite entanglement or violation of Bell inequalities. The requested 
measurement fidelity is already available for charge qubits, where significant progress 
has been achieved with dispersive readout inside a cavity, providing a visibility of more 
than 90% [ID]. A similar design has been proposed for flux qubits [H]. Moreover, other 
experiments based on Josephson junction technology indicating similar fidelities have 
been performed [321 031 01 051 E] ■ 



Efficient creation of multipartite entanglement in flux qubits 13 




1 0.95 0.9 0.85 0.8 0.75 0.7 

Fidelity f 

Figure 7: Expectation values of optimal Bell operators for GHZ type, W type, and 
CHSH type entanglement shown as a function of decreasing measurement fidelity. 

Table 2: Minimal detector fidelities for the detection of tripartite entanglement or 
violation of Bell inequalities, respectively. 



Operator f min Operator /, 
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6. Related work 

Related work has shown similar properties for a ring of exchange-coupled qubits [46] 
even in the ground state. Open/linear coupling topologies, albeit easier to prepare 
experimentally, require more complex pulse sequences [47j|48j|49] because the eigenstates 
do not have tripartite entanglement; they become more efficient in connected networks 
[50] . Also, tripartite entanglement between two superconducting cavities and one qubit 
has been proposed [HI]. Beyond tripartite entanglement, a circuit QED setup has been 
suggested for the fast preparation of an iV-qubit GHZ state in superconducting flux or 
charge qubits [52] . 
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7. Conclusions 



In this paper, we introduced a system of three coupled flux qubits in a loop. We showed 
that it exhibits strong tripartite entanglement for a realistic and approachable set of 
parameters and that it is possible to detect and quantify this entanglement. 
We acknowldege support by the Deutsche Forschungsgemeinschaft through SFB 631 and 
the NSERC disovery grants program and, in parts, by the National Science Foundation 
under Grant No. NSF PHY05-51164. 



Appendix A. State preparation in a degenerate subspace 

We describe here, how the highly entangled states (e.g., I^max)) m the degenerate 
subspaces introduced in section [3] can be prepared through driving by a resonant laser 
field. Let \g, n) denote the ground state of the system dressed by |n) photons and 
1^1/2, n — 1) the excited states with one photon less. This constitutes a ^-level scheme. 
Under resonant driving, these states are all degenerate but can be coupled using the 
effective transition Rabi frequencies u\/2 (which depend on the coupling to the field) as 



to 



red 



h 








w 2 e 







-up 



\ 



) 



where the phase ip can be introduced by a relative phase between the driving fields. 
Diagonalization of this Hamiltonian leads to eigenstates that are hybridized between the 
coupled qubit system and the driving photons. One finds that under the corresponding 
time-evolution operator the application of a pulse of length t = 7r/2ya>f + u\ to the 
ground state \g, n) = (1, 0, 0) leads to the final state 



i 



i ,2 

cu 2 



\ 



Using multiple coils, arbitrary combinations of u\ and u 2 as well as relative phases can 
be produced and thus, arbitrary linear combinations of the degenerate states can be 
prepared. This type of preparation requires multiple, phase-locked microwave drives, 
ideally originating from the same source. This type of phase stability is routinely 
achieved in high-precision microwave control [53l 154"] . As this is essentially a precise 
Rabi pulse, it is relatively robust against decoherence. However, after the preparation 
of the state, decoherence in the degenerate subspace essentially selects the eigenstates 
of the coupling to the environment, projected to the subspace, which will in general be 
not or much more weakly entangled. Hence, although the preparation of these states is 
as robust as any preparation of eigenstates, their maintenance is, in general, as fragile 
as that of a superposition. 
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Appendix B. Local Decomposition of entanglement witnesses 

Table IB1I lists the decompositions of the used entanglement witnesses into local 
projective measurements. The decompositions are optimal in that they require the 
minimal set of local measurements (e.g. four measurement settings for VV^^: af 3 -i.e. 
a z on all qubits-, a® 3 , (a z + a x )® 3 , and (a z - a x )® 3 ). 

Table Bl: Local decomposition of the entanglement witnesses used above. The 



GHI 

)® (four settings) was calculated from = «x(^)® 3 W^ 

GHZ GHZ GHi- 

decomposition for WW (five settings) was obtained similarly from the optimized 
decomposition W^ (^ ve settings) derived in [33], WW = U L ^ VV^ U L ■ 



optimal decomposition for can be found in (four settings), the one for 

w£L (four settings) was calculated from W^L = R x (tt)® 3 W^LiU-vr)® 3 . The 

GHZ GHZ GHZ 



EW Local decomposition 



ng^ ± [10 • H® 3 + 4 af - 2(a y a y t + a y ta y + t a y a y ) - (a z + <r x ,)® 3 - (a z - a x )® 3 ] 
i [10 • II® 3 - 4 of 3 - 2(a y a y t + a y la y + l a y a y ) + {a z - a x ) m + (a z + o^)® 3 ] 



WW ± [17 ■ ll® 3 - 7cr® 3 - 3(6", 11+ ta x 1 + tta x ) + 5(a x a x 1 + a x 1 a x + 1 a x a x ) 
-(1 - a x + a z ) <g> (1 - a x - ^ a y - \ a z ) <g> (1 - a x + ^ a y -\b z )- 
-{t-a x - a z ) <g> (11 - a x + ^ a y + \ & z ) ®{t-a x -^§ a y + \& g )- 
-(11 - b x + a y ) <g> (1 - b x - \ a y + ^ <r z ) ® (11 - - i o-j, - ^ o- z )- 
-(11 -a x -a y )®(t-a x + l a y - ^ £ z ) <g> (1 - a x + \ a y + ^ <r 2 )] 
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